Prom the Hubbard to the SO (5) Ladder: A Numerical Study 
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The Hubbard Hamiltonian on a two-leg ladder is studied numerically using quantum Monte Carlo 
and Exact Diagonalization techniques. A rung interaction, V , is turned on such that the resulting 
model has an exact 50(5) symmetry when V = —U. The evolution of the low energy excitation 
spectrum is presented from the pure Hubbard ladder to the SO(5) ladder. It is shown that the low 
energy excitations in the pure Hubbard ladder have an approximate SO (5) symmetry. 
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In recent years, various compounds containing two-leg 
cuprate ladders have been synthesized Q]. The ground 
states of these materials are characterized by strong elec- 
tronic interactions with a variety of competing ground 
states. Theoretically, the focus has been on using the 
Hubbard [Q and t — J models || to understand the 
ground state properties of these materials. At half-filling, 
the two-leg ladders exhibit a spin-gap insulating phase, 
while away from half- filling d x 2_ y 2 pairing and charge 
density wave (CDW) correlations become enhanced 0. 

Recently, Zhang proposed a model with SO (5) symme- 
try to relate the two phases of antiferromagnetism (AF) 
and <i-wave superconductivity (dSC), even though these 
two phases have quite distinct ground states ||. Sub- 
sequently, several groups have constructed microscopic 
models which have an exact 50(5) symmetry Fur- 
thermore, it has been argued that both the Hubbard |l(]] 
and t — J models jll| in two-dimensions (2D) have ap- 
proximate 50(5) symmetry. 

In this paper, we examine one such microscopic Hamil- 
tonian H using quantum Monte Carlo (QMC) and exact 
diagonalization (ED) techniques. The model we study 
can be changed adiabatically from the pure Hubbard 
model on a two-leg ladder to a model with 5*0(5) symme- 
try, thus elucidating how the low energy spectrum evolves 
as the higher symmetry is approached. In this way, we 
can test whether the low lying excitations of the Hubbard 
ladder are well described by the 5*0(5) picture. 

The microscopic 50(5) ladder we investigate here was 
first introduced by Scalapino, Zhang and Hankc (SZH) 
||. The SZH model contains only local interactions on a 
rung of the ladder and is therefore much easier to visu- 
alize and implement numerically than 2D 50(5) models 
with long range interactions. The Hamiltonian is given 
by 
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Hhop allows electrons to move along the ladder (in the x 
direction) as well as within a rung (in the y direction): 
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where c] creates an electron at site i with spin projection 
s. The interaction term, Hi nt , is given by 

flint = Uj^inii - l/2)(n u - l/2) + ^n is 

i i,s 

+y$>i-l)(n i+ y-l) + J^5V5Vy (3) 

i i 

where n\ is the number operator for site i; U is the 
Coulomb repulsion between electrons occupying the same 
site; V is an interaction between electrons on a given 
rung; J is a rung spin-spin interaction with 
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In order for this Hamiltonian to be manifestly 50(5) 
symmetric, the condition of J = A(U + V) must be im- 
posed. 

The phase diagram of the SZH model at half-filling has 
been obtained at strong coupling, i.e., U, V >> t±,tn 
B. Predominately for U > 0, the ground state is well 
described as a product of rung singlets, i.e., a spin-gap 
insulator. Doped holes in this ground state form d x 2_ y 2 
rung pairs and exhibit power law pairing and CDW cor- 
relations. A phase transition occurs for positive values 
of U from the singlet ground state into a triplet ground 
state along the line V = —U, where each rung is occupied 
by a triplet magnon or a doublet pair. Furthermore, the 
triplet ground state can be considered to be an 50(5) 
generalization of the spin-one Heisenberg chain, i.e., a 
system with a finite excitation gap and short range cor- 
relations. 

At weak couplings (U,V << tx,tn), the phase dia- 
gram of the SZH model was obtained using a perturba- 
tive renormalization group analysis ]l2[ ]. The resulting 
phase diagram, although consistent with the strong cou- 
pling result, contains new phases not found in the strong 
coupling limit. Furthermore, it was shown that two-leg 
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Hubbard-like ladders flow to 50(5) symmetry even when 
explicit symmetry breaking terms were included, such as 
longer range hoppings O], 

In this paper, we explicitly turn off the spin-spin inter- 
action by requiring that J = E3]. Thus, our model is 
50(5) symmetric only when V = —U, where the ground 
state at large couplings was found to have a degeneracy 
between the rung singlet and triplet 50(5) multiplets. 
Consequently, we now have an adjustable model which 
interpolates from the Hubbard ladder to the 50(5) lad- 
der by varying V from to —U. We will study this 
model in order to investigate how the low lying excita- 
tions evolve as V — > — U, i.e., as the model approaches 
50(5) symmetry. Throughout, we will take the isotropic 
hopping case of t = t± = tu = 1 and will only consider 
U > 0. 

In Zhang's 50(5) theory, the superspin vector takes 
on a fixed magnitude below some characteristic temper- 
ature, T*, and its direction fluctuates between the AF 
and dSC phases. For the case in which the model is 
manifestly SO (5) symmetric (V — —U), T* is pushed to 
infinity and all of the components of the superspin vector 
will be equal. To test this, we use deterministic quantum 
Monte Carlo (QMC) Jljj] to measure the first and fourth 
components of the superspin vector given by 
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{~ic\r)a y c f (r + y) + h.c) (5) 



and 



n 4 (r) 



(-1)' 



(c\r)o- z c(r) - c f (r + y)cr z c(r + y)). (6) 



Here r is the rung index, and the spinor index on the c 
operators has been suppressed, i.e., c(r) = (c| )r , cj,, r )- 
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FIG. 1. (ni) (circles) and (714) (squares) for U = At and 
f3t = 2 at half-filling for the 2x8 (open symbols) and the 
2 x 12 (filled symbols) lattice. 

In Fig. 1, the sum of ri\ and over all rungs of 
a 2 x 8 and 2 x 12 ladder are shown, i.e., («i(4)} = 



jv^£>I(4)W"i(4)(0)>- At half-filling (/1 = 0) for 
an intermediate strength coupling of U = At, the ex- 
pected behavior is seen: as — V approaches the SO (5) 
value of U, the correlations become equal. Note that 
even at half-filling a sign problem exists due to the new 
Hubbard- Stratonovich fields introduced by the V term 
p6| . At small to intermediate coupling strengths, i.e., 
U < 6t, reliable results can be obtained, even at the 
50(5) point when V — —U. However, when U becomes 
large, and consequently —V becomes large as the sym- 
metric point is approached, the sign problem becomes 
unmanageable. Thus, only U = At with (3t = 2 QMC 
results are shown |l7| . 

Recall that at strong couplings along the phase tran- 
sition line of U = —V, the ground state was found 
to have a degeneracy between the rung singlet and 
rung triplet 50(5) multiplets. Therefore, the corre- 
lation functions should decay as a power law in dis- 
tance, r, with some thermal exponential activation, i.e., 
n i(4)( r 7/3) ~ exp(— (3/{;)/r a (where £ is the thermal cor- 
relation length). However, when V is different from —U, 
the correlations should decay exponentially in r. 
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FIG. 2. (a) m(r) = (n|(r)m(0)} for U = it and fit = 2 at 
half-filling on a 2 x 12 ladder as a function of distance along 
the ladder for three different values of the rung attraction V . 
Note that the results for the SO(5) ladder are for V = —4t 
(open diamonds), (b) Same as in (a) for ni{r) = (n\(r)ri4(p)) . 

In order to test if this critical behavior could be seen, 
the correlations of the superspin components were mea- 
sured as a function of distance along the ladder. In Fig. 
2, we show the results obtained for the 2 x 12 ladder 
with U — At at half-filling as V approaches the 50(5) 
value. In Fig. 2(a), a qualitative change in the correla- 
tions of ni(r) is clearly observed between the Hubbard 
ladder and the 50(5) ladder, indicating a crossover from 
a power law to an exponential decay in n\{r) as the 50(5) 
point is approached. In contrast, the change in the n^ir) 
correlations (Fig. 2(b)) is less pronounced. 

It was not possible to cleanly extract the power law 
exponents from a log-log analysis of these correlations 
for several reasons. First, U = At is not strong enough to 
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recover the exact degeneracy between SO(5) multiplets 
which is present only at large coupling strengths. Also, 
at a relatively high temperature of fit = 2, the thermal 
correlation length tends to dominate the behavior of the 
superspin components j2j. Since reliable measurements 
at either larger couplings or lower temperatures could not 
be obtained, let us now turn to a zero-temperature exact 
diagonalization (ED) analysis to better understand the 
low energy spectrum ]fq] . 



What is the nature of the ground state of the 2x6 
ladder? In Fig. 3(a), the dynamical spin response, using 
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FIG. 3. (a) Exact diagonalization results for the dynamical 
spin response for a 2 x 6 ladder at half-filling with U = At as 
V goes from zero up to the SO(5) value of —At. (b) Same 
as (a) for U = 8t and V varying from zero to — St. (c) The 
spin-gap as a function of U for the SO(5) case with V = —U. 
The solid line is the strong coupling result, J = At 2 /U. 

Because of the large Hilbert space at half-filling, the 
system sizes which can be studied with ED are limited, 
and we will only present results for a 2 x 6 ladder. On the 
other hand, an advantage of ED is that it allows for the 
direct calculation of the spin-gap, as well as dynamical 
spectra 10 nj given by: 



0(q,w) 



-Im(k|OL, 



7T 



lu + Eg, s , — H — it 



Oq|k) 



(7) 



where |k) is the ground state wave function with mo- 
mentum k, and 0(q, uj) is an arbitrary translationally 
invariant operator. 
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S* = J2j exp(— iq • r^S^rj), is shown for U 
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|]l9| . Starting from the Hubbard ladder, V — 0, we find 
that a spin-gap is present at half-filling 0. In particu- 
lar, as seen in Fig. 3, as V approaches — Z7, the spin-gap 
remains robust and even increases for the SO (5) ladder 
rather than going to zero as one might expect from the 
strong coupling (J — > 0) limit. The low lying spin excita- 
tions are dominated by a large peak in the spin response 
occurring at a momentum transfer of (tt, tt). (Note that 
Figs. 3(a) and (b) show the momentum integrated spin 
excitation spectrum, S(u>) — J dqS'(q, uj).) Furthermore, 
there is little qualitative difference between the spectra 
as the model moves toward the ladder of higher symme- 
try. Similar behavior is seen for U — 8t (Fig. 3(b)) where 
again a spin-gap is always present. An angular resolved 
examination of 5*(q, uj) shows that the dominant low en- 
ergy peak in S(uj) is due to inter-band scattering pro- 
cesses, involving a momentum transfer of tt in the rung 
direction, whereas the quasi-continuous spectral weight 
at higher frequencies stems from intra-band processes. 

The spin-gap can also be obtained by measuring the 
energy difference between the ground state in the S z — 
sector and the ground state with total S z = 1 at a mo- 
mentum of {tt, tt), shown in Fig. 3(c) with V = —V. 
The value obtained in this way agrees exactly with the 
position of the lowest peak in S(uj). The gap increases, 
reaching its maximum around the intermediate coupling 
strength of U — 6t, and then falls off as 1/U. For com- 
parison, the expected result from strong coupling given 
by At 2 /U is shown in Fig. 3(c) as the solid line. There- 
fore, it is understandable why the power law behavior in 
the correlations was not seen from the QMC simulations. 
An effective spin-spin interaction, J e ff = At 2 /U, causes 
the system to have a spin gap with the ground state be- 
ing well described by a product of rung singlets. Hence, 
no degeneracy between the two SO(5) multiplets exists. 

In order to explore the approximate symmetry of the 
Hubbard ladder and how it evolves as a function of V, 
the dynamic response of the tt operator, i.e., O q = 7r q = 
^ r exp(— iq ■ r)-7r(r) with 

Tr(r) = i(^(r)+^(r)) = ^(-l) r c T (r + y)c T (r), (8) 

was measured for the state with one electron pair more 
than the half-filled state. Fig. 4 shows a direct compar- 
ison at U — At and U = 8t between the low energy spin 
excitations and the tt resonance using the energy of the 
half-filled ground state as a reference. Clearly, both have 
a dominant low energy peak at the same excitation en- 
ergy independent of the rung coupling V, indicating that 
the resulting final states are equivalent. Thus, an approx- 
imate SO(5) symmetry is revealed. However, it should 
be noted that the tt resonance has significant spectral 
weight at higher energies (not shown in Fig. 4) for the 
Hubbard ladder, whereas all the weight shifts to the low 
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energy peak as the symmetric ladder is approached p0| . 
A thorough finite-size scaling analysis is expected to show 
that this low energy ir resonance survives in the thermo- 
dynamic limit as the length of the ladder is increased. 
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FIG. 4. (a) Comparison between the low lying spin exci- 
tations of a 2 x 6 ladder at half-filling with U — At and the 
7r excitation spectrum of the ground state with one pair of 
electrons above half-filling, (b) Same as (a) for U = 8t. Note 
that the energies are taken with respect to the ground state 
energy at half-filling. 

In this paper, we have numerically analyzed a model 
which can be varied by changing a rung interaction so 
that it interpolates between the Hubbard and the 5(3(5) 
SZH model on a two-leg ladder. We find that the spin-gap 
vanishes as U ft goes to infinity, and the ground state is 
well described by rung singlets for intermediate to large 
coupling strengths. This can be understood since the 
hopping along the chain and within a rung creates an 
effective spin-spin interaction J e // = At 2 /U. This effec- 
tive spin-spin interaction must be overcome by an even 
stronger negative value of the rung interaction, —\V\ > U 
in order to restore the degeneracy between rung singlets 
and triplets obtained as U/t — > oo. Furthermore, no 
qualitative changes occur in the low lying excitations of 
the Hubbard ladder as the 50(5) symmetric point is ap- 
proached. Therefore, we conclude that the 7r operators 
are approximate eigenoperators of the Hubbard ladder 
near half-filling, and hence, the low energy behavior of 
the Hubbard ladder in the intermediate to strong cou- 
pling range has an approximate 50(5) symmetry. 
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